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We describe a versatile toolbox for the quantum simulation of many-body lattice models, capable of exploring 
the combined effects of background Abelian and non-Abelian gauge fields, bond and site disorder, and strong 
on-site interactions. We show how to control the quantum dynamics of particles trapped in lattice potentials 
by the photon-assisted tunneling induced by periodic drivings. This scheme is general enough to be applied to 
either bosons or fermions with the additional advantage of being non-perturbative. It finds an ideal application 
in microfabricated ion trap arrays, where the quantized vibrational modes of the ions can be described by a 
quantum lattice model. We present a detailed theoretical proposal for a quantum simulator in that experimental 
setup, and show that it is possible to explore phases of matter that range from the fractional quantum Hall effect, 
to exotic strongly-correlated glasses, or flux-lattice models decorated with arbitrary patterns of localized defects. 
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I. INTRODUCTION 

In general, quantum many-body systems cannot be under- 
stood by extrapolating the properties of their individual con- 
stituents [1], but rather by considering their rich collective 
behavior. A prototypical example is that of sound waves in 
solids [2], where the atoms do not move individually but vi- 
brate collectively and give rise to propagating quasiparticles, 
the so-called phonons. This particular many-body problem is 



exceptional since the properties of the phonons can be calcu- 
lated exactly, something that rarely occurs in quantum many- 
body systems. The usual paradigm is that the complexity 
of the model grows very fast with the number of particles, 
making both analytical and numerical methods more involved 
and less efficient. A radically different approach is based on 
the so-called experimental quantum simulations, which were 
envisaged several decades ago by R. P. Feynman [3], and 
have now evolved into a discipline that merges concepts from 
atomic physics, quantum optics, quantum-information science 
and condensed-matter physics. The central idea of a quantum 
simulation is to manipulate the microscopic properties of a 
particular experimental setup in a way that it reproduces faith- 
fully a quantum many-body model under study. In this way, 
nature itself computes the properties of the model, and our 
measurements yield the answer to questions such as the na- 
ture of the ground- state and collective excitations, determin- 
ing whether a dedicated model is sufficient to describe the rel- 
evant properties of a particular system. Accordingly, quantum 
simulations have the potential of solving fundamental open 
problems in physics, ranging from high-temperature super- 
conductivity, to the thermalization of closed quantum systems, 
or the properties of spin glasses. 

Imagine for a moment that the phonons in a solid were con- 
trollable to such an extent that their Hamiltonian could be 
engineered to target a many-body model of interest. For in- 
stance, by shaping the crystal anharmonicities, one could tune 
the phonon-phonon interactions and reproduce the physics of 
strongly-correlated models. Unfortunately, it is hard to de- 
velop techniques that allow such microscopic control in solid- 
state materials, and it might be beneficial to search for dif- 
ferent experimental platforms. This exotic idea may become 
realized in experiments with Coulomb crystals of cold atomic 
ions in radio-frequency traps [4]. Current experimental tools 
allow for a promising control of the vibrational excitations at 
the quantum level, such that phonon-based quantum simula- 
tions of many -body physics may become a reality [5, 6]. 

So far, the most exploited property of the vibrational exci- 
tations of trapped ions is their ability to transfer information 
between distant ions, which encode quantum bits (qubits) in 
their electronic states. This provides a mechanism to perform 
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quantum logic operations between distant qubits, a funda- 
mental building block of quantum-information protocols [7] . 
This type of two-qubit couplings can be interpreted as a spin- 
spin interaction where the qubit plays the role of a pseu- 
dospin (henceforth referred as spin) [9] and can be exploited 
for quantum-simulation purposes. From this perspective, the 
phonons act as mediators of the interaction and give rise to 
a wide range of spin models [10], which are familiar in the 
field of quantum magnetism. The experimental success of 
these spin-based quantum simulations, either in an analog [11] 
or digital version [12], has motivated a variety of propos- 
als that range from neural networks [13], to three-body in- 
teractions [14], frustrated magnetism [15], mixed-spin mod- 
els [16], or topologically-ordered spin models [17]. 

Instead of using the phonons as a gadget to obtain the de- 
sired quantum-spin simulator, one can substantially enhance 
the capabilities by reclaiming phonons as the building blocks 
for the quantum simulation of many -body models [5]. The 
local vibrational excitations of each atomic ion, considered to 
be trapped individually, correspond to bosonic quasiparticles, 
and the Coulomb interaction is responsible for the interchange 
of these bosons between distant ions. To build interesting 
bosonic quantum simulators, one must complement the afore- 
mentioned scheme with additional ingredients, such as strong 
trapping non-linearities that lead to Mott insulating [5] and 
frustrated phases [18], or incommensurate trapping potentials 
leading to different quantum phase transitions [19]. Besides, 
the spins of the ions can be used as a tool to widen the appli- 
cability of this quantum simulator, which can potentially tar- 
get the spin-boson model [20], the lattice Jay nes -Gumming s 
model [21], or the phenomenon of Anderson localization due 
to disorder [22]. In addition to these many-body quantum sim- 
ulations, there has also been a recent activity in building ana- 
logues of single-particle phenomena with a special emphasis 
on relativistic effects (see e.g. [23]). In this case, even if the 
models are tractable on a classical computer, the predicted ef- 
fects are hard to access in the original experiments, and thus 
justify the effort in building a trapped-ion analogue. 

In this article, we introduce a versatile toolbox for the quan- 
tum simulation of many-body lattice models. This toolbox is 
based on the phenomenon of photon-assisted tunneling (PAT) 
of phonons in ion traps [24], but can also be applied to dif- 
ferent systems of bosons, and even fermions, in a lattice. As 
shown in this manuscript, the PAT effect has a a variety of 
facets that can be exploited to provide new paradigms for the 
aforementioned many -body quantum simulation. The idea un- 
derlying the PAT is that the tunneling of particles between the 
wells of a periodic lattice can be assisted by inducing reso- 
nances that correspond to the absorption/emission of photons 
out of an EM-field providing a periodic driving force. Even 
though this idea was initially introduced for condensed-matter 
systems [25] (see [26] and references therein), it has also been 
applied in the field of ultracold atoms in optical lattices [27] . 
Here, the analogy is straightforward since the neutral atoms 
can tunnel between the adjacent wells of a driven periodic po- 
tential created by light. In this work, we show how to exploit 
the PAT in order to induce synthetic gauge fields in the afore- 
mentioned many-body lattice models. 



In the context of trapped ions, the possibility of control- 
ling the tunneling of the vibrational excitations between ions 
stored in two separate potential minima, aligned within a lin- 
ear radio-frequency trap, has been recently demonstrated [28]. 
The exchange of phonons between two ions can be con- 
trolled and optimized by matching their individual trapping 
frequencies. Building on these seminal experiments, we have 
demonstrated theoretically that driving the trapping frequen- 
cies in/out of resonance realizes a novel version of the PAT 
paradigm [24]. Rather than relying on a periodic force [25, 
27] scaled to a larger number of ions, which may be quite de- 
manding in trapped-ion experiments, our scheme relies on the 
periodic modulation of the trapping frequencies. In particular, 
it requires an additional relative phase between the individual 
modulations that can be experimentally tuned. This allows 
for the full control of both the amplitude and the phase of the 
tunneling of phonons, a result that becomes specially interest- 
ing for two-dimensional arrays of micro-fabricated traps [29]. 
Here, the technology of ion-trap micro-fabrication provides 
a means of assembling arrays of ion traps in any desired ge- 
ometry [30]. The combination of the capabilities for the de- 
sign and fabrication of arrays of microtraps with the tool of 
PAT [24] opens a vast amount of possibilities for a quan- 
tum simulator of many-body physics. These range from the 
realization of bosonic models subjected to Abelian and non- 
Ahelian synthetic gauge fields, to bond and site disorder lead- 
ing to glassy phases, or to decorated flux lattices which can 
be related to anyonic excitations. In this article, we provide a 
detailed analysis of such PAT in arrays of microtraps. 

This article is organized as follows. The main results are 
summarized in Sec. IV, so that this section can be consulted 
for a general overview. The general scheme of PAT for any 
lattice model is introduced in Sec. II, where we demonstrate 
how the synthetic gauge field arises from the assisted tunnel- 
ing. The application as a toolbox for quantum simulations 
with trapped ions is contained in Sec. III. Here, we present a 
thorough description of the different many -body models that 
can be explored by exploiting the peculiarities of trapped ions, 
which widen the versatility of our quantum simulator. 



II. PHOTON-ASSISTED TUNNELING AND SYNTHETIC 
GAUGE FIELDS 

In this Section, we describe the concept of PAT for a general 
system, and discuss how it can be exploited for the quantum 
simulation of lattice models incorporating synthetic gauge 
fields. Once the main ingredients are identified, we focus 
on the case of ions in micro-fabricated traps in Sec. Ill, and 
show how to implement the PAT for the vibrational excita- 
tions. Let us remark that the generality of this section may 
also find applications in different platforms, such as photons 
in arrays of cavities in circuit QED [31], or ultracold atoms in 
optical lattices [32]. We also note that the PAT scheme is not 
restricted to bosons, but works equally well for fermions. This 
may alleviate some of the difficulties that arise in the simula- 
tion of synthetic gauge fields via Raman-assisted tunneling of 
fermionic atoms in spin-dependent optical lattices [33]. 
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Figure 1: Photon-assisted tunneling scheme: (a) Scheme for the 
driven one-dimensional tight-binding model in Eqs. (1) and (4). The 
tunneling of the particles between neighboring sites r? ^ rP, which 
is initially suppressed by the large gradient Aft)cj ^ becomes as- 
sisted by a resonant periodic modulation of the on-site energies with 
amplitude y]^(0^. (b) Two-dimensional scheme giving rise to a syn- 
thetic gauge field such that the path around an elementary plaquette 
oc q^^o mimics the Aharonov-Bohm phase, which is picked up by 
a charged particle following the path around a plaquette pierced by 
an external magnetic field orthogonal to the lattice. 



We consider a tight-binding model describing the tunneling 
of particles, either bosons or fermions, between the sites of an 
underlying two-dimensional lattice. The lattice sites are char- 
acterized by a vector of integer numbers i = (11^2), such that 
r? = /iJiCi -\- 12(^2^2, where {ta} are the unit vectors span- 
ning the lattice, and {da} are the corresponding lattice con- 
stants. The particles are represented by bosonic or fermionic 
creation-annihilation operators j where a labels some 

additional degrees of freedom. The dynamics of the system is 
described by the following Hamiltonian 

H = Ho^H, = Y, COaA^^a^ + L L ('^taj4,i^cTj + H.c.) , 

CT,i o i>j 

(1) 

where cOfj i stands for the on-site energy (h=l), and J^^^ is the 
tunneling amplitude of the particles between different lattice 
sites j ^ i, which usually depends on the distance between 
sites such that /^^j = /^^dr? — r?|). The tunneling constraint 
i > j refers to an ordering of the sites, such that ii > j\ , or 
h > h if h = ji- Two additional ingredients are required: 

i) Gradient of the on-site energies: The on-site energies 
have the following expression cOa i = COa -\-AcO(jii. Here, COa 
is a constant energy offset, and AcOa results due to a gradient 
along one of the lattice principal axes satisfying 



(2) 



ii) Periodic modulation of the on-site energies: The on-site 
energies must be supplemented by the periodic modulation 
(^a,i (^a,i + Vd(^d cos{(0^t + 0i), where (Od is the driving fre- 
quency, and r](i(Od the driving amplitude. Note that the peri- 
odic driving incorporates a site-dependent phase 0i that shall 
play a crucial role for the PAT, and is described by 



These two ingredients, schematically represented in 
Fig. 1(a), are incorporated in the above description by modi- 
fying the Hamiltonian Hq Ho{t) as follows 

^o(0 =^(cOcy+A(Ocy/i +r]da)dCos((Od^ + 0i))a^ia^i. (4) 



Let us emphasize again that the standard formulation of the 
PAT relies on a periodic force acting on the particles residing 
on the lattice sites [25, 27]. In contrast, our scheme consid- 
ers a periodic modulation of the on-site energies, which turns 
out to be better suited to taylor the tunneling amplitudes in 
analogy with a background synthetic gauge field. In subsec- 
tion II A, we present an analytical model, which is confronted 
to numerical simulations in II B. 



A. Analytical description 

1. Dressed tunneling 

In this subsection, we derive a compact expression for 
the PAT strength, and study its dependence on the periodic- 
driving parameters. Let us express the tunneling Hamiltonian 
Hi in the interaction picture with respect to Ho{t), namely 
Ht{t) = U{tyHtU{t) where U{t) = Q-ifo^^^oiT) ^ ^^lis pic- 
ture, the annihilation operators fulfill 



Acta A 



((Ocj+A(Ocj/i +r]d(OdCos(a)d^ + 0i))aa,i, (5) 



which leads to the following relation between both pictures 



^(7,1- 



(6) 



At this point, we note that the tunneling Hamiltonian Hi is in- 
variant under U(l) gauge transformations, and thus, the last 
term in the above expression can be trivially gauged away 
^cT,i e~^^dsin((^i)^^ . Accordingly, the Hamiltonian becomes 

Htit) = Ii>j-^d';ij(Oi,i«<.j +H.C., where J^.-it) = J^^^it), 
and the time-dependence is encoded in the function 



0(f) 



,iT7d(sin(a)d?+0i)-sin(a)d?+0j)) iAcOa(ii-ii> 



Jr7d (sin(a)d?+(^i)-sin(cOd?+0j)) 



if ^'i > 71 
if ii 



J\- 
(7) 

In this expression, one readily observes that the tunneling am- 
plitude becomes dressed by the photons of the periodic driv- 
ing J^^^ -^^ijlO- Besides, the fundamental role of the phase 
of the periodic driving also becomes apparent: only when 
0i 7^ 0j, the tunneling becomes assisted /^^ (t) ^ J^^y In order 
to proceed with this analytical treatment, we use the identity 



^iT7dsin(a)d?+0i) _ 



(8) 



0i = 01/1+02^2. 



(3) 



where Jsi^ld) are Bessel functions of the first class [34]. 
Hence, the dressed tunneling becomes a sum of terms with 
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different time dependences. Let us first focus on the tunneUng 
along the gradient, ii > j\ , which can be expressed as follows 

(9) 

By tuning the driving frequency to (0^ = ^(Oq / r, where r is 
some positive integer, the above expression contains two dif- 
ferent types of terms. There are resonant terms fulfilling 

/ = ^ + r(/i-7i), (10) 

whereas the remaining terms are far off-resonance. By apply- 
ing a rotating- wave approximation (RWA) for J^^^ <C AcOfy (2), 
we neglect the rapidly-oscillating terms yielding 

Jh = Jml.J^W,^ri,H-n){r],^^^^^^^ (11) 

According to this expression, both the amplitude and the 
phase of the tunneling are controlled by the periodic driving 
parameters. As shown in Sec. Ill, the periodic driving can be 
achieved by means of optical dipole forces on the ions, and 
the tunneling along the ei-axis is thus assisted by the absorp- 
tion/emission of r(i\ — ji) photons from the EM-field provid- 
ing the driving force, hence the name PAT. A similar analysis 
for the tunneling orthogonal to the gradient, /i = ji, shows 
that the resonant terms satisfy = s. 

The complete assisted-tunneling Hamiltonian becomes 

^eff = LL^'^ij4,i^CTj +H.C., (12) 

(J i>j 

with the dressed couplings expressed as follows 

•/d'ij =-/5j^/(iJ)(^d,r/d,A<^.ij)e-'^(*+^j), (13) 

where the function ^ is responsible for the modulation of the 
tunneling amplitude 

<^z(C,^,0) = l^MQJs+xi^y^'^^^', (14) 

seZ 

and we have defined the phase difference A^y = 0; — 0j , and 
the function /(ij) = r{h - ji). 

Let us emphasize one of the important properties of this 
photon-assisted scheme, namely, its non-perturbative charac- 
ter. So far, we have not used any restriction on the values of 
r]d or 0i, so that the dressed tunneling may be of the same 
order as the bare one J[. In Fig. 2, we represent the modulation 
function ^/(ij) (T7d, T7d, A0ij) for different periodic-driving pa- 
rameters, and different ranges of the tunneling. In particular. 
Figs. 2(a) -(c) represent the tunneling along the direction of 
the gradient for r = 1, whereas Fig. 2(d) stands for the tun- 
neling orthogonal to the gradient. As derived from Fig. 2(a), 
for T7d ~ 1 , A0 = TT, one gets the maximal amplitude for the 
assisted tunneling between nearest neighbors. In particular, 
we find 1 ^ 0.6 for these optimal values, which thus sup- 
ports our previous claim about the non-perturbative nature of 
the scheme (i.e. \J^.^j^^^ J ^ 0.6|/^^j^^^ J). It is also interesting 




Figure 2: Photon-assisted modulation of the tunneling ampli- 
tude: Contour plot of the modulation amplitude for the 
tunneling between sites i j, such that ii > ji, as a function of 
the driving parameters T]^, A^y for r = 1. (a) First-neighbor assisted 
tunneling (i ^ i + ei), (b) Second-neighbor assisted tunneling (i 
i + 2ei), (c) Third-neighbor assisted tunneling (i ^ i + 3ei), (d) First 
neighbor assisted tunneling orthogonal to the gradient i ^ i + e2. 

to note that, as we consider longer range terms, their ampli- 
tude gets gradually diminished, as shown in Figs. 2(b)-(c). In 
Fig. 2(d), we represent the modulation amplitude in the di- 
rection perpendicular to the gradient, which taking the same 
parameters as above leads to \ J^.[^q^ J ~ 0.2|/^*^j^^^ J. 

Before concluding this subsection, we briefly comment on 
another interesting property of this scheme. Let us consider 
that the original Hamiltonian (1) also contains on-site particle- 
particle interactions H-\-V, where 

gg' i 

and the dependence of the interaction strengths U(j(j' on the in- 
ternal indices depends on the bosonic/fermionic nature of the 
particles. It is straightforward to see that the gradient and pe- 
riodic driving do not modify V. Therefore, the PAT also works 
in the presence of on-site interactions. It would be interesting 
to study how the scheme can be used to modify long-range in- 
teractions, but this lies beyond the scope of the present work. 

2. Synthetic gauge fields 

In this subsection, we demonstrate that it is possible to in- 
terpret the phase of the dressed tunneling (13) as if it was orig- 
inated by a background gauge field. In particular, we show 
that whenever particles tunnel along a closed path in the lat- 
tice, they pick up a non- vanishing phase analogous to the cele- 
brated Aharonov-Bohm phase for charged particles in electro- 
magnetic fields [35]. The consecutive tunneling of a particle 
around a unit plaquette of the lattice i^i + ei^i + ei+e2^ 
i + 62 ^ i [Fig. 1(b)] is formally given by 

11/(1) _ 

^^O ~ '^d;i,i+e2'^d;i+e2,i+ei+e2'^d;i+ei+e2,i+ei'^d;i+ei,i- ^^^^ 



Using Eqs. (7) and (13), it can be expressed as 



(1) 



|/t^(J2)/t''(^i)^o(r]d,T]d,02)^r(r]d,r]d,0i) 

(17) 

which leads to W^^ = |WQ^^|e^^^, where 0o = r02 only de- 
pends on the component of the periodic-driving phase that is 
orthogonal to the gradient. This quantity is a gauge-invariant 
observable proportional to the so-called Wilson loop in lattice 
gauge theories [36], which can be expressed as follows 



(1) 



•dS 



(18) 



Here, we have introduced an effective charge ^* independent 
of the charged/neutral character of the particles, the synthetic 
gauge potential As, and the synthetic gauge field Bs. Inde- 
pendently on how we rearrange the phases of the tunneling 
strengths locally, the enclosed phase 0o = ^* /□ Bs • dS shall 
always be left invariant, and can be thus interpreted as the 
magnetic flux of a synthetic magnetic field Bs that pierces the 
lattice. In order for this analogy to be complete, one should 
consider carefully the long-range character of the tunneling, 
according to which particles can follow different closed paths. 
We shall explore the two possible paths around the second 
smallest plaquettes, namely 



(2) 



/(2) 



r<7 jG jO jO 

' '^d;i,i+2e2'^d;i+2e2,i+ei+2e2'^d;i+ei+2e2,i+ei'^d;i+ei,i' 

r<7 jG jG jG 

' '^d;i,i+e2'^d;i+e2,i+2ei+e2'^d;i+2ei+e2,i+2ei'^d;i+2ei,i- 



(19) 



Repeating the above calculations, we find that /|^( 



(2)/lw(2)|_ 



^O^Vl^O^^I = Hence, the Aharonov-Bohm phase is 

doubled with respect to the unit plaquette, which is consistent 
with the fact that the enclosed area is also doubled for these 
paths. The same occurs for any other closed path, and thus 
the analogy to a background gauge field is consistent with the 
possible long-range of the tunnelings. 

At this point, it is worth emphasizing the generality of 
the scheme hereby proposed. It works for both bosons and 
fermions, for any range of the tunneling J^^^ = /^^(|r? — r?|), 
it can incorporate local particle-particle interactions, and it is 
non-perturbative. In the second part of this manuscript, we 
shall specify this scheme to trapped-ion experiments, which 
provide an ideal platform where to realize a bosonic PAT (i.e. 
vibrational excitations). Before moving onto the numerical 
verification of these results, let us introduce an alternative and 
more compact formulation of the synthetic gauge fields. It 
is possible to rewrite the PAT Hamiltonian in analogy to the 
so-called Peierls substitution [37], which yields 



G i>j 



(20) 
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Figure 3: Photon-assisted tunneling along a link: (a) Bare tun- 
neling for an initial state |^o) = ^^jlvac) with a single particle occu- 
pying the first site. The numerical expectation values {ni{t)) (blue 
circle), and (^2(0) (orange squares) are compared to the effective 
description (wi(O)eff (blue lines), and (w2(0)eff (orange lines) de- 
scribed in Eq. (21). (b) Suppressed tunneling due to the gradient, 
and in the absence of periodic driving, (c) Assisted tunneling in 
the presence of a periodic driving for A0i2 = 7i. (d) Coherent de- 
struction of tunneling (black arrows) due to a periodic driving in the 
absence of the gradient. The yellow circles represent the maximum 
population transferred to site 2 = max{(w2(0) • ^ < ^*}^ such that 
t* = lOOn /\^(){ri(^,7r)\co, as obtained from the numerical integra- 
tion of Eqs. (1) and (4). The red line represents the predictions for 
the same magnitude ^2 according to the effective Hamiltonian (12). 



Xi = 01 ^1 Z-^- ^^^^ ^^^^ ^hi^ transformation is also consistent 
with the long-range character of the tunneling. This compact 
formulation (20) condenses the main result of this section: a 
gradient and a periodic modulation of the on-site energies give 
us full access to the amplitude and phase of the tunneling, 
such that the effects of background synthetic gauge fields can 
be mimicked even for neutral particles or quasiparticles. 



B. Numerical simulations 

In this subsection, we confront the analytical model for the 
PAT in Eq. (20) with numerical simulations for the tunneling 
Hamiltonian (1) subjected to the gradient and periodic driv- 
ing (4). We analyze in detail some basic realizations of the 
PAT, which allow us to carry out a thorough numerical study 
considering even the effects of finite temperatures. 



where we have introduced the dressed-tunneling amplitude 
^^ij = -^tUj^/liJ) (^d, ^d, A0ij), and the synthetic gauge poten- 
tial As(x) = —Boyei, where = r02/^*<^i<^2- This gauge po- 
tential, which corresponds to the famous Landau gauge for 
electrons in a constant magnetic field, is obtained after the 
following U(l) gauge transformation aa^i aa^iQ~^^\ where 



1. Photon-assisted tunneling along one link 

Let us initially focus on the simplest situation to test the ac- 
curacy of the analytical model, namely, a lattice consisting of 
only two sites populated by spinless bosons. We consider the 
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following parameters of the driven Hamiltonian in Eqs. (1) 
and (4): the bare tunneling is /t;i2 = 10~^ CO, the gradient 
Aft) = 0.5 CO, and the periodic driving parameters correspond to 
r]d = 1, and co^ = Aco. Note that we do not consider additional 
degrees of freedom, avoiding thus the index a in the follow- 
ing. In Fig. 3, we represent the expectation values {ni{t)) = 
that result from the numerical integration of 
the Schrodinger equation id\'¥{t))/dt = {Ho{t) +ift)|^(0)' 
and compare them with the effective analytical description 
1^(0 )eff = e~'^eff^|^o). We consider an initial state with a 
single bosonic particle in the first site |^o) = ^^jlvac), where 
I vac) stands for the vacuum. In this case, the effective descrip- 
tion (12) can be integrated exactly, yielding the following pe- 
riodic oscillations of the particle between the two lattice sites 

(^l(O)eff = ^(l+COS(OeffO, ^^^^ 
(^2(0)eff = ^(l-COS(OeffO, 

where (Ogff = 2|/d;i2| is the frequency of oscillations. In the 
absence of any gradient or periodic driving, cOgff = 2|/t;i2| = 
0.02a), and thus the particle exchange undergoes oscillations 
with a period of Tt-n = ^/\Jv,i2\ = lOOn/ CO, which coincides 
exactly with the scale shown in Fig. 3(a). In this figure, we re- 
veal that the effective description (circles and squares) agrees 
with the numerical results (solid lines). Note that in order to 
carry out the numerical integration, we truncate the bosonic 
Hilbert space to ^tmnc = 4. We have confirmed that the trun- 
cation error after changing ^tmnc '^tmnc = ^tmnc + 1 lies be- 
low |(n/(0)^trunc - (^K0)/2truncl < l^"! Thc samc occurs in 
Fig. 3(b), where we switch on the gradient. As a consequence 
of the mismatch between the on-site energies, the tunneling 
is completely inhibited. In order to assist it, we switch on 
the periodic driving in Fig. 3(c) for 0i = 02 + ^, which acti- 
vates again the periodic oscillations but with a different period 
Tt-n = ^/\Jd;i2\ = 1007r/|^i(r7d,r7d,7r)|a). For the parame- 
ters used, we obtain |^i(T]d,T]d,^)| ^ 0.6 (Fig. 2(a)), which 
explains the slightly longer oscillation period. However, we 
remark that both dynamics occur on the same time-scale as a 
consequence of the non-perturbative character of the scheme. 

In Fig. 3(d), we study the maximal population transfer to 
site 2 due to the periodic driving r]d ^ 0, but in the absence 
of the gradient Aco = 0. We set the phases to 0i = 0, 02 = 
and study the maximal transfer = max{(^2(0) • ^ ^ ^*}' 
such that ^* = 100n/\^o{ri^,ri^,n)\co for a range of driv- 
ing strengths r]d G [0,4]. As shown in this figure, there are 
certain values of the driving strength, marked by black ar- 
rows, where the tunneling gets completely suppressed, a phe- 
nomenon known as coherent destruction of tunneling [25], or 
dynamic localization. Let us remark that this scheme leads 
to a perfect localization of the particles also in a longer one- 
dimensional chain, since the dressed tunneling cancels for all 
pairs of nearest-neighboring lattice sites simultaneously. Once 
we set = ni\/2, such that A^^-^^^+i = tt, it is possible to find 
a zero of the modulating function ^o(^d ' ' ~ ^ which 
is related to the particular values of the Bessel functions. In 
this limit, the particles shall not diffuse through the lattice. 

In Fig. 4(a), we study the photon-assisted process for dif- 
ferent driving phases A0i2 G [0, In], considering the same pa- 



rameters and initial state as before. We represent the max- 
imum population transferred to site 2 ^2 = max{(n2(0) • 
t < f], such that f = 50n /\^i{ri^,ri^,7r)\o) corresponds 
to the optimal PAT with A0i2 = 71. The yellow dots corre- 
spond to the numerical integration of the truncated Hamil- 
tonian (1) and (4), whereas the solid line corresponds to the 
analytical prediction evaluated at the related exchange pe- 
riod ^2 = (^2(^0)eff for the different phases, namely = 
SOtt/I^i (T7d, T7d, A012) |(0. In this figure, we reveal the agree- 
ment between both descriptions for any of the driving phases. 
Note that this agreement will not be compromised by going to 
larger arrays. As announced previously, we can interpolate be- 
tween the optimally-assisted and the completely-suppressed 
tunneling regimes by modifying the driving phase. 

It is also interesting to consider an initial state that does not 
correspond to a single particle, but rather to a thermal ensem- 
ble. We focus now on bosons with independent mean number 
of particles /zi , /22- This state corresponds to po = Pi P2 with 

p. = 1 ^ = ^^(^^-^ima]at^^ ^ ^22) 

such that the parameters ft are implicitly defined through the 
mean number of particles following a Bose-Einstein distribu- 
tion hi = l/(e^'^' ~ !)• Ill this case, we must integrate nu- 
merically the Liouville-Von Neumann equation idp{t)/dt = 
[Ho{t) -\- Ht^p{t)], and compare the result to the effective de- 
scription Peff(0 = e~'^eff^pQe+^^eff^ In particular, we con- 
sider the same parameters as before, and set fii = 0.5, /22 = 
0.25 after truncating the particle Hilbert space to ^tmnc = 7. 
Note that due to the thermal effects, the truncation parame- 
ter has been increased with respect to the previous simula- 
tions. Again, we have checked the convergence of the re- 
sults by increasing ^tmnc- In Fig. 4(b), we represent the ex- 
pectation values ^2 5 ^2 introduced above, together with those 
of site 1, nl = min{(^i(^)) : r < ^*}, and n[ = (^i(^O)eff for 
t' = 5O;r/|^i(T7d,T7d,A0i2)|co. We conclude from the perfect 
agreement shown in the figure that the scheme also works for 
thermal states. In fact, in the absence of interactions, the equa- 
tions of motion do not depend on whether the state is pure or a 
mixed thermal state. For the latter, there is a background over 
which the PAT phenomena will occur. The role of the addi- 
tional interactions, and its interplay with the thermal states, is 
an extremely interesting topic that deserves a separate study. 

The suitability of the PAT scheme for thermal states will 
turn out to be essential in Sec. Ill, where we discuss a realis- 
tic implementation of the quantum simulator of gauge fields 
with phonons in microtrap arrays. In this case, it means that 
cooling to the vibrational ground-state is not necessary to ob- 
serve the non-trivial effects of the PAT. Another interesting 
topic is the occurrence of motional heating in the microtraps. 
In case this motional heating is global, it should not inter- 
fere with the PAT effects. However, if this heating has a local 
nature, it may lead to very interesting effects that mimic the 
role played by a reservoir providing electrons to local regions 
of a metallic conductor. This effect can be complemented by 
the controlled engineering of dephasing in the phonon dynam- 
ics by introducing noisy potentials in the trap electrodes [38]. 
Hence, the applicability of the PAT quantum simulator could 
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Figure 4: Phase-dependence and thermal photon-assisted tun- 
neling: (a) Dependence of the PAT on the relative driving phase 
A012. The yellow circles represent the maximum population transfer 
to site 2, fZ2, for a particle initially located at site 1, |^o) = |vac), 
as obtained from the numerical integration of the full time-dependent 
Hamiltonian HQ{t) +//t. The red line corresponds to the same quan- 
tity evaluated for the effective description //gff' referred as Em- 
phasized by the doted frame is the region of optimal assisted tunnel- 
ing, (b) Phase-dependence of the PAT for an initial thermal state with 
fii = 0.5, = 0.25. We also represent the particle transfer to site 1. 

be widened, and used to study how transport phenomena is 
affected by an environment that induces decoherence effects 
such as motional heating or dephasing. 



2. Photon-assisted tunneling around a plaquette 

Let us now apply our formalism to a two-dimensional setup 
where to test the application of PAT for the quantum simula- 
tion of synthetic gauge fields. We consider a model of spin- 
less bosons in a square plaquette, and set the parameters of 
the total Hamiltonian in Eqs. (1) and (4) to /t;i2 = 10~^a), 
Aco = 0.5(0, r]d = 1, r = 1, and ft)d = Aco. Our aim is to test 
the accuracy of the effective description (20), trying to find 
clear-cut evidence of the underlying synthetic gauge field. In 
order to do so, we study a discrete version of the celebrated 
Aharonov-Bohm effect [35]. 

Let us briefly recall the underlying interference effect 
[Fig. 5(a)]. A particle initially localized at site 1 can tunnel 
to site 3 following two different paths, namely 7i : 1 ^ 2 ^ 3, 
or 72 : 1 ^ 4 ^ 3. Due to the gauge field, each path leads to a 
different phase |^f) oc e''* k ^^'^^a^lvac) +e''*^r2d''-^^a^|vac), 
such that the probability of performing such trajectory is 
P3 - 2 + 2cos (^* /^^ dr • As - ^* 4 dr • As) = 2(1 + cos 0o), 
where 0o = ^* dr • As was defined in Eq. (18). From this ex- 
pression, one finds a perfect destructive interference that for- 
bids the particle to tunnel to site 3 when the enclosed phase 
00 = TT. This Aharonov-Bohm interference will serve us as a 
testbed for the validity of the gauge-field analogy (20). 

We now compare the dynamics obtained from the numer- 
ical integration of the driven Hamiltonian (1) and (4), to the 
effective description in Eq. (20). In Figs. 5(b)-(c), we con- 
sider l^o) = ^^jlvac), and study the propagation of the par- 
ticle along the square plaquette for different driving phases. 
In Fig. 5(b), we observe that for the phases 0i = TT, (^2 = 0, 
the particle is allowed to tunnel to every lattice site. This is 



consistent with the fact that for 0^ = 02 = 0, there is no in- 
terference. Conversely, we set 0i = ;r, 02 = ^ in Fig. 5(c), 
where one readily observes that the tunneling to site 3 is for- 
bidden by the aforementioned Aharonov-Bohm interference. 
We note that for the trapped-ion case (see the sections below), 
the diagonal path going from sites 1^3 and 2^4 must also 
be accounted since the tunneling amplitude need not be small. 
In [24], we discussed how to cancel it by playing with the 
phases, so that the analogy with the Aharonov-Bohm effect is 
valid. In Fig. 5(e), we corroborate that the effect also holds for 
thermal states. Finally, we check in Fig. 5(d) that the perfect 
destructive interference only occurs for the so-called TT-flux 
phase. In this figure, we represent the maximal transferred 
population to site 3, ^3 = max{(^3(r),0 < ^ < r*}, where 
^* = 1007r/|^i(r7d,r7d,7r)|a), and show that the perfect inter- 
ference only occurs when the enclosed phase 0o = ^- This 
allows us to rule out other possible sources of interference, 
and conclude that it is only due to the synthetic gauge field. 

Let us close this subsection by highlighting the accuracy of 
our analytical treatment, which has been confronted to a nu- 
merical survey of different setups with a wide range of param- 
eters. The rest of this manuscript will be devoted to analyze 
the experimental setups in the field of quantum optics, ions 
in microtrap arrays, where the necessary ingredients for this 
PAT are in reach with state-of-the-art technologies. Besides, 
we shall also introduce novel methods to to exploit additional 
aspects of PAT phenomenology. 



III. PHOTON-ASSISTED-TUNNELING TOOLBOX FOR 
TRAPPED IONS 



In this section, we apply the above scheme to arrays of 
micro-fabricated ion traps, and show that one can achieve PAT 
of the vibrational excitations between distant microtraps. In 
order to set the notation, we start by showing in Sec. Ill A 
that the vibrations of an ensemble of ions in a microtrap ar- 
ray can be described in terms of tunneling phonons [5]. Then, 
we describe in Sec. IIIB how the gradient and periodic driv- 
ing of the trapping frequencies can be achieved using the 
tools of state-of-the-art experiments in ion traps (see the re- 
views [46]). In particular, the gradient can be implemented 
by the local control of DC-voltages applied to the trap elec- 
trodes, whereas the periodic driving stems from an optical 
dipole force. We derive a set of constraints that this dipole 
force must fulfill, and discuss how these requirements can be 
met for realistic experimental parameters. In Sec. Ill C, we 
incorporate phonon-phonon interactions in the phonon-based 
quantum simulator of bosonic particles exposed to synthetic 
gauge fields. Finally, in Sees. HID, HIE, and IIIF, we de- 
scribe how to extend the PAT scheme beyond the applications 
discussed so far, reaching more involved models where the 
spin of the ion can be exploited as an additional tool. 
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Figure 5: Photon-assisted tunneling around a plaquette: (a) 

Schematic representation of the plaquette, highhghting the two pos- 
sible paths connecting sites 1 and 3. (b) PAT for an initial state 
l^o) = <3||vac), such that the periodic -driving phase is 0i = tt, and 
0(3 = r02 = with r = \. The numerical expectation values {ni{t)) 
(blue circles), (^2(0) (^^^ diamonds), (n^it)) (green stars), and 
{n4{t)) (yellow squares) are compared to the effective description 
in Eq. (20) (solid lines with same colors). For vanishing 0q = 0, 
the particle can tunnel anywhere in the lattice, (c) Same as above, 
but setting (pQ = tt, where one clearly observes the Aharonov-Bohm 
interference (^3(0) = 0- (d) Maximum population transfer to site 3 
as a function of the enclosed phase 0o E [0,27r]. The yellow dots 
correspond to the exact numerical integration n^, and the solid red 
line to the effective description ^3 . A black arrow marks the perfect 
Aharonov-Bohm interference, (e) Aharonov-Bohm interference for 
an initial thermal state fii =0.25, and n2 = = n4 = 0.1 



A. Tunneling of phonons in microtrap arrays 

Let us consider an ensemble of N ions with charge e and 
mass m, labelled by integer numbers i = (ixjy)- Each ion is 
subjected to the electric potential energy created by a planar 
micro-fabricated trap Vt, and the Coulomb energy potential Vc 
causing the repulsion between the remaining ions, such that 
the total Hamiltonian is 



(23) 



Here, the trap is designed in such a way that the ion equilib- 
rium positions, which are given by V(yt + Vc) |j,o = 0, form a 

regular lattice r? = /iJiCi -\-i2d2e2, where ei,e2 are the unit 



vectors, and ^1,^2 the lattice spacings. At sufficiently small 
temperatures, the excursions of the ions from these equilib- 
rium positions rj = rP + Srj are small with respect to the lattice 
spacing. Thus, for our purpose, we are allowed to consider the 
above Hamiltonian up to second order only [2, 47]. This ap- 
proximation yields a system of coupled harmonic oscillators 



2m 



-EI(C + ^c-/)5'-«,i5r,j, 



(24) 



with the following couplings 
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where rP_j = r? — 5[\ stands for the Kronecker delta, and 
a, 7 = X, refer to the main axes of the trapping potential. 
Note that the micro-fabricated trap gives rise to a confinement 
that can be considered to be harmonic for the motional am- 
plitudes considered here, and is therefore characterized by the 
frequencies (Oa i that may depend on the axis and the ion posi- 
tion within the lattice. The Hamiltonian (24) can be expressed 
in the basis of local quantized vibrations 



^'"«''=V2i::K'+^0' 



(26) 



where b^a\,b^^ are the bosonic creation-annihilation opera- 
tors, and we set /i = 1. In this local basis, the vibrational 
Hamiltonian becomes H = Hq + Hf^, 

a,i ij «,r 

(27) 

where the Coulomb couplings /^^ = / {2m-^(0a^i(0Y^^) de- 
scribe the exchange of vibrational excitations between differ- 
ent ions i,a ^ j,7, and yield the aforementioned collective 
phonons once diagonalized. Let us express in the interac- 
tion picture with respect to Hq, namely H'^ = e^^^^ HcQ~^^^\ 



i{c0a^i-C0y^j)t ^ ^t^ ./?t ^e''(f«a,i+^7j)^)+ H.C. 



In this work, we consider that the microtraps fulfill 



(28) 



CO, 



I ? 

(29) 

which allow us to neglect the rapidly oscillating terms of the 
above Hamiltonian using a standard rotating-wave approxi- 
mation (RWA). The first condition allows us to consider inde- 
pendent vibrations along each of the trapping axes, whereas 
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the second one allows us to neglect the terms that do not con- 
serve the number of vibrational excitations. Accordingly, 

H = Ho + H, = + 11 +H.C.), 

i,a oc i>j 

(30) 

where the trapping frequency of each ion is slightly modi- 
fied by the electrostatic interaction with its neighboring ions 
<,i = (^l,i^^cT/^^y^^^ and we have defined /^^jj =2/^^^. 
This Hamiltonian can be interpreted as a tight-binding model 
for the local phonons, which tunnel between distant micro- 
traps according to the dipolar tunneling strengths /^j. By a di- 
rect comparison to the general tight-binding model in Eq. (1), 
we can identify the index a with the vibrational axis, the tun- 
neling strength J^^^ with the aforementioned dipolar Coulomb 
couplings, and the on-site energies cOfj i with the microtrap 
trapping frequencies. In the following subsection, we describe 
how to obtain the gradient and the periodic driving using state- 
of-the-art tools in trapped-ion experiments. 



B. Gradient and periodic driving of the trapping frequencies 

In microtrap arrays, it is possible to design the individual 
trapping frequencies {cOa.i} by the control of DC-voltages ap- 
plied to local micro-fabricated electrodes. In fact, the experi- 
ments [28] have made explicit use of this property in a conven- 
tionally segmented linear rf-trap, in order to switch on/off the 
phonon tunneling by tuning the trapping frequencies of two 
neighboring traps in/out of common resonance [29] . We chose 
the trapping frequencies distributed according to a gradient 
(^a,i = (Oa-\- ^COah' Considering mutual trap and mutual ion 
distances > 10/im due to current constraints in fabrication. 
The correction to the trapping frequency due to the electro- 
static interaction can be neglected, and we get the desired gra- 
dient of the on-site energies 0)^ j ^ COa i = C0a-\-A(0aii which 
must be incorporated to the tight-binding Hamiltonian (30). 

We now turn to the periodic driving of the trapping fre- 
quencies, which shall be implemented by an optical dipole 
force [46] . We consider a pair of laser beams with frequencies 

, ft>2, and wavevectors ki , k2, which couple to the electronic 
states I ii), I ti)^ l^i) (see the energy-level scheme of Fig. 6 for 
the particular case of ^^Mg+ ions). Encoding the spin de- 
gree of freedom into a pair of electronic ground states by ex- 
ploiting their hyperfine or Zeeman splitting, the energy differ- 
ence (Oq/Ik between the states | ii) , | ti) lies in the microwave 
range 1 GHz, whereas the energy splitting (O^/ln to an auxil- 
iary state |ai) is much larger (Oa ^ (t^, and typically lies in the 
optical domain 100 — 1000 THz. By tuning the beatnote of 
the phase-locked laser beams to a particular value, it is possi- 
ble to choose out of a variety of couplings. For instance, when 
(0\ — ayi^(OQ, such that coi , 0)2 <C (Da, one speaks about a stim- 
ulated two-photon Raman transition between the ground and 
excited states [Fig. 6(b)]. Conversely, when a)i — 0)2 ~ (O^ <C 
(Oq, one obtains a running- wave realization of the so-called 
spin-dependent dipole forces. In this work, we are interested 
in yet a different regime (0\— ayi<^(Oa [Fig. 6(c)], where we 
get the desired periodic driving by the laser dipole force. We 
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Figure 6: Energy-level scheme and laser-beam arrangement: (a) 

Reduced diagram for the hyperfine energy levels |F,M) of'^^Mg^. 
The ground-state manifold is split into the hyperfine levels F = 
2,3, such that an additional magnetic field allows us to isolate two 
magnetic sublevels | ti) = |2,2)i, | |i) = |3,3)i with an energy differ- 
ence of ft)o/27r ^ 1.8 GHz. Note also that these levels are dressed by 
a vibrational ladder with equidistant spacing (Oa/^Ti ^ 1 — lOMHz. 
The excited manifold /2 contains an auxiliary level to implement 
a stimulated Raman scheme with two laser beams and a beatnote 
(0]^ = (Oi — 0}2. (b) Raman scheme for two laser beams with Rabi fre- 
quencies Cli , Q.2 in the so-called red-sideband regime CO]^'^ — (Oa- 
In this regime, excitations are exchanged between the electronic and 
vibrational states of the ions, (c) Raman scheme for the periodic 
modulation of the trapping frequencies COl <C CD^. In this limit, only 
negligible coupling between the internal or vibrational states take 
place during the transition, but rather a periodic modulation of the 
trapping frequencies is realized. 



consider the ion-laser interaction in the dipolar approximation 
^L= 1^ ^^^li?|ai)(ti|e-('^^«)+H.c. 

/=1,2 i 

XXi^l«|ai)ai|e-('^^r,--^^^ 



(31) 



/=1,2 



where Q.^^}. , Q.^^ are the Rabi frequencies for each of the transi- 
tions to the auxiliary level induced by each of the laser beams. 
Since we are assuming that the laser frequencies are far off- 
resonant with respect to the auxiliary state, it is possible to 
perform an adiabatic elimination of this state to simplify the 
dynamics between | Furthermore, for sufficiently 

small Rabi frequencies and considering that (0\ — Oi <^ (Oq, 
it is possible to obtain 



i 

-I 



l^Vl^gi(Ak-ri-a)LO +H.C. 



4A 



4A 



gi(Ak-ri-a)LO ^H.C. 



i;i)aii 



iti)(tii, 



(32) 



where we have defined the beatnote (Ol = a)i — 0)2, Ak = 
ki - k2, and A = (Oa - ftio - a)i such that |A| > (Oq. Note 
that we have neglected the AC- Stark shifts arising from each 
single laser beam, since they do not play any role for the 
periodic driving of the trapping frequencies. By a proper 
choice of the laser intensities, detunings, and polarizations, 
one finds a regime where the two-photon Rabi frequency is 
= -(a[\^Ya^^^/2A = -(a[i^Ya^^^ /2A, and thus the 



^L = -KT^i:72A=-K^ 

laser-ion interaction becomes 



Hl = iS^LliC'^^'''"'-'*^') +H.C. 



(33) 
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Table I: Typical parameters for ion microtrap arrays. 
1-10 MHz 50-500 kHz 50-500 kHz 1-iokHz 0.1-1 MHz 0.1-0.4 



This expression corresponds to a Stark shift that acts equally 
on both electronic states, and is caused by the crossed laser 
beams of different frequencies. 

To obtain the desired periodic driving from the ion-laser 
Hamiltonian (33), we express the ion position in terms of the 
local phonon operators ri = rP + Y.a (^^ i + i) / V^2ma)^. 
By introducing the so-called Lamb-Dicke parameter r{a = 

• Ak/ y/2m(0a <C 1, we can Taylor expand the Hamiltonian 

i^L ~ ^L^^ ^ +^L^^ whereby 

^0)^^^^^,Ak.rO-.-a^.^H.C. 
^ i 

'my 

+ H.C. 

(34) 

Note that the first term corresponds to an irrelevant c-number 
that does not modify the phonon dynamics. In order to find the 
relevant contribution of the two remaining terms H^\h^\ 
we switch to the interaction picture with respect to Hq = 
-fia i' we consider that both, the frequency gra- 
dient and the laser frequency, fulfill AcOf^ , (Ol <C (Oa^\(Oa — 
(!^j\ai-y^ a rotating wave approximation for r[a^\. <^ al- 
low us to neglect the rapidly-oscillating terms and leads to the 
desired periodic driving of the trapping frequencies 

i^L - -^lLl^r7^cos(Ak-rP-a)L0/^^Vi- ^^^^ 

By direct comparison with Eq. (4), we find the following, in- 
dividually tunable, mapping between the parameters of the pe- 
riodic driving and the laser beams: cOd = (Ol, (OdT]^ = — ^^L^a, 
and 0i = — Ak • rP. Let us also note that in order to neglect the 
higher order terms that rotate with the laser frequency (Ol, we 
impose that ^^l^^ It is important to remark that even if 

we are considering small Lamb-Dicke parameters r7of <C 1, the 
periodic driving |r]d| = ^^L^^/^i^ needs not to be small by se- 
lecting the appropriate laser and Rabi frequencies, so that the 
non-perturbative character of the PAT can be exploited. 

The remaining task to demonstrate that the scheme of PAT 
works for phonons in microtrap arrays is to consider the 
typical values for the experimental parameters, and discuss 
whether they satisfy the constraints imposed during the above 
derivation. The orders of magnitude of all the relevant pa- 
rameters are listed in Table I, which satisfy the different con- 
straints made along this derivation, namely 

r7a^lL<C0a, ACO^, (Ol <C (Oa, |(0a - (Oy|a^y. (36) 



For the last inequality, it suffices to focus on the phonon 
modes transverse to the microtrap plane, a = z, such that 
Ico^ — ft)y|y^^/27r ^ IMHz. If the in-plane vibrational modes 
are to be used, the anisotropy of the trapping frequencies 
(Ox ^ (Oy should be considered in the microtrap design. Fi- 
nally, table I allows us to estimate the parameters of the model 
Hamiltonian, and reveal whether the constraints for the PAT 
in Eqs. (2) and (3) are fulfilled. We find that J^^^l^^ ^ 1" 
10kHz<C AcOa/27: ^ 50-500kHz, which thus satisfies the con- 
straint in Eq. (2). Besides, we find that 0i = 0i/i + 02^2, where 
0a = — (Ak • ea)da, and thus the constraint (3) is also satis- 
fied. Therefore, we can conclude that the required ingredients 
for the PAT introduced in Sec. II can already be met with the 
current technology of microtrap arrays if the tunneling rates 
outrun the decoherence rates. Note that if motional heating 
turns out to degrade the results severely, we could mitigate 
this problem by the using a cryogenic setup. 

Let us finally comment on two additional constraints differ- 
ent from (2) and (3), which are particular to the trapped-ion 
setup. Both the gradient and the the periodic driving should 
not modify the stability of the ion crystal, and thus must be 
smaller than the trapping frequencies 

AcOa,ri^co^<^cOa. (37) 

This is also fulfilled for the parameters in Table I. Let us note 
that this condition sets a limit to the scalability of our pro- 
posal. For gradients A(Oa/27l ^ 50kHz, and trapping frequen- 
cies cOa/^TT ^IMHz, an array of 10x10 microtraps is still 
consistent with the maximum attainable trapping frequency. 
We note that the gradient can be reduced further (lOkHz) with- 
out compromising the efficiency of the PAT scheme, and lead- 
ing to larger systems with N =2500 microtraps. Finally, we 
would like to remark that the scheme could be scaled even 
further by considering a local gradient that only affects a few 
sites, and repeats periodically along one axes of the micro- 
trap array. In order to assist the tunneling between the sites 
where the gradient is changed, while maintaining a homoge- 
nous synthetic flux, an additional periodic driving with a suit- 
able frequency and phase must also be introduced. 

C. Synthetic gauge fields and phonon-phonon interactions 

The discussion of the PAT of phonons has focused on 
the quadratic Hamiltonians corresponding to the Coulomb- 
induced tunneling (30), and the periodic driving of the trap 
frequencies (35). However, as discussed in Sec. II, non- 
linearities corresponding to the on-site interactions in Eq. (15) 
can be incorporated without modifying the assisted- tunneling 
scheme. In fact, this broadens the applicability of our many- 
body quantum simulator, since it also targets models of 
strongly-interacting particles. To obtain such phonon-phonon 
interactions [5], strong non-linearities in the trapping potential 
are required. Remarkably, one can exploit a different realiza- 
tion of the above dipole forces (33) in order to give rise to such 
non-linearities. We denote the parameters of these new dipole 
forces by a wiggled bar in order to differentiate them form the 
dipole forces leading to the periodic driving. 
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Table II: Mapping of microtrap phonons onto the PAT model 
Eqs. (30),(38): b^. a /«y (Ol -^liI Ak-rf ^^Cli^fflfi^cos^i 



We consider that two additional laser beams leading to a 
dipole force have the same frequency, in a way that the beat- 
note d)L = 0, and Eq. (33) corresponds to a standing wave. 
In analogy to the previous expansion (34), we consider that 
^hVa (^a so that all the terms that do not conserve the 
number of phonons can be neglected. However, since d)L = 0, 
the quadratic terms now correspond to a small shift of the trap- 
ping frequencies which has no effect since ^^l^^ AcOa <C 
COa. The most relevant contribution is due to the quartic terms 

i,a,7 

(38) 

which corresponds exactly to the on-site interactions intro- 
duced in Eq. (15). In Table II, we summarize the mapping of 
the phonon Hamiltonian onto the original periodically driven 
tunneling Hamiltonian of Sec. II. As discussed there, since 
this interaction is purely local, it shall not be modified by the 
scheme of PAT. Therefore, the effective phonon Hamiltonian 
becomes H^ff = K^ff + Vgff , where 

" (39) 

where the tunneling is = ^Hj^/(i,j)(^d7 ^d^^^ij). and the 
gauge field can be expressed as As(x) = —Boyei, such that 
Bq = (r02/^*Ji J2). Equation (39) incorporates the central re- 
sult of this section, which shows that the PAT in microtrap 
arrays leads to a quantum simulator of a long-range Bose- 
Hubbard model [42] under additional gauge fields [43]. 

There are several interesting regimes for this quantum sim- 
ulator. In the non-interacting limit for a square microtrap ar- 
ray, it yields a bosonic dipolar version of the so-called Azbel- 
Harper-Hofstadter model [39]. The nearest-neighbor model 
has been studied thoroughly during the last decades, and con- 
tains several interesting properties that range from its relation 
to topological numbers, to the fractal and self-similar proper- 
ties of its energy spectrum and wavefunctions, the existence of 
gapless edge excitations, or the so-called TT-flux phases [40]. 
The addition of the long-range dipolar tunnelings introduces 
a new feature in the model that, to the best of our knowledge, 
has not been studied previously and may modify the above 
phenomena. Besides, most of these effects rely on a magnetic 
flux per plaquette on the order of the flux quantum, which 
cannot be achieved in solid-state materials assuming realistic 
magnetic fields. On the contrary, our proposal has the poten- 
tial to reach these regimes since it is non-perturbative and the 
flux can attain arbitrary values 00 G [0,2;r]. 

Since it is possible to build any desired microtrap geom- 
etry [30], our quantum simulator can explore the physics of 



bosonic ladders subjected to synthetic gauge fluxes. Besides, 
the capability of tuning the fluxes, together with the indepen- 
dent control of the tunneling strength along/across the ladder 
rungs, dives into the phenomena of flat-band physics and edge 
states [41], which is typical for fermionic topological insula- 
tors that break the time-reversal symmetry [45]. 

Another interesting regime corresponds to ^ Ui^ay^ 
where the interactions compete with the kinetic energy and 
induce strong correlations in the Bose-Hubbard model [42]. 
With respect to the neutral-atom realizations [44], the phonon 
model includes the effects of longer-range tunnelings, and the 
possibility to address site-dependent interactions. Besides, 
thanks to the tunability of the functions ^/(ij) [Fig. 2], the 
strongly-correlated regime can be reached, in principle, re- 
gardless of how small the on-site interactions are. Let us 
note, however, that there is a fundamental limit to this ap- 
proach, which is imposed by external sources of decoherence 
and heating. Accordingly, the dynamics must always be faster 
than the time-scale imposed by these sources of noise. 

Eor vibrational modes transverse to the microtrap array 
a = z, the condition Sl\^f\a <C (Oa can be relaxed for ions in the 
node of the standing wave. It then suffices to set Sl\^f\a <^ 
which allows us to reach interaction strengths in the ~1-10 
kHz-regime, which are directly on the order of the bare tun- 
nelings. Finally, including the gauge fields in this strongly- 
correlated regime further enhances the versatility of our quan- 
tum simulator. Even if the particles are bosonic and the in- 
teractions are local, one can target fractional quantum Hall 
states, and composite-fermion fluids [43]. 



D. Non-Abelian synthetic gauge fields 

The synthetic gauge fields discussed so far (39) correspond 
to standard electromagnetism. In this theory, they are for- 
mally introduced to restore the invariance with respect to local 
unitary transformations in the gauge group U(l). A natural 
question that arises is whether these synthetic fields can be 
generalized to different gauge groups, possibly non-Abelian 
ones [48]. In this case, the local unitary also acts on some ad- 
ditional degree of freedom, which we shall refer to as the fla- 
vor. Here, we show that it is possible to realize such scenarios 
by exploiting two or three orthogonal directions of vibration 
as the different flavors of the non-Abelian theory. 

We first exploit two main axes of vibration a=x^y within 
the plane defined by the microtrap array, the corresponding 
trapping frequencies being different (0^^ (Oy. In the regime 
where the bare tunneling strength is smaller than the fre- 
quency difference, these two directions are uncoupled (see 
Eq. (30)). The main idea is to use two gradients of opposite 
sign for each degree of freedom, and two separate periodic 
drivings of the same frequency but of different amplitudes 

(^x,i = (0x^ Acoil + ri^co^ cos{co^t + 0i) , ^^^^ 
cOyi = cOy-Acoii + 77dya)dcos((0d^ + 0i). 

According to Eq. (35), this modulation of the trap frequencies 
can be obtained from a single Raman-beam scheme, such that 
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the corresponding wavevector has a component along both 
axes r7djc = —^lVx^ Vdy = —^hVy- repeating the anal- 
ysis of Sec. II, and considering the same type of conditions, 
we reveal that the dressed tunneling strengths (13) must be 
modified for each of the vibrational axes 

(41) 

where /a(i,j) = ra{i\ — Ji), and one must introduce the fol- 
lowing axis-dependent parameter = r^Vy = —r. It is then 
straightforward to rewrite the effective Hamiltonian accord- 
ing to a generalized Peierls substitution 

i>j 

(42) 

where the synthetic gauge potential now depends on the re- 
lated vibrational axis, and is equivalent to A^(r) = — As (r) = 
—Boy^i, such that = r02/^*<^i<^2- Let us introduce now a 
bosonic spinor field operator = to describe the 

phonon fields corresponding to vibrations in each direction. 
The kinetic part can be written 

/^eff = lyiK^''^^'^^^ +H.C., (43) 

i>j 

where K^-^ is a matrix that describes the vibrational couplings 
for each spinor component, 

and A^^ = —B^yXz^i is a SU(2) non-Abelian gauge field, that 
we write in terms of a Pauli matrix, T^, acting on the vibra- 
tional index (i.e. flavor index). The associated magnetic field 
corresponds to B^^ = BoT^e^, so that each flavor is subjected 
to an opposite flux piercing the lattice. For fermions, these 
types of gauge fields give rise to the so-called quantum spin 
Hall effect [49], which is the prototype of time-reversal pre- 
serving topological insulators in two-dimensions [45]. Let us 
note that Eq. (43), together with the arguments presented in 
Eq. (19) imply that each flavor around a plaquette accumu- 
lates a non-Abelian Aharanov-Bohm phase that is governed 
solely by the SU(2) gauge field A^^. However, in addition to 
that, each flavors, y, is subjected to different tunnelings, /jj j, 
^.j j , such that an spin-orbit coupling is superimposed to the 
non-Abelian gauge. The latter effect may enrich the dynamics 
with respect to the usual situation in SU(2) gauge theories. 

We remark that the above synthetic gauge field (43) is only 
a particular type of non-Abelian SU(2) gauge fields. In order 
to consider more general fields, it is possible to exploit the 
quadratic terms of the Coulomb interaction (28) that mix the 
vibrational modes along x^y. By setting the frequency of the 
periodic driving to account for both the frequency difference 
(cOx ^ COy) and the particular gradient, the dressed tunnelings 
would also involve a change of the flavor index, yileding thus 
more general gauge fields in the group SU(2). 



E. Spin-mediated disordered Hamiltonians 

The properties of solids usually differ from those of per- 
fectly periodic crystals. In realistic samples, there is a certain 
amount of disorder in the form of impurities, dislocations, or 
vacancies, which may alter dramatically the properties of the 
solid. The study of disorder in solids is an active and ma- 
ture field of condensed matter [50], where the system Hamil- 
tonians are usually modeled as stochastic operators. Here, 
the randomness is due to a statistical description of the dis- 
ordered degrees of freedom. To incorporate such a random- 
ness in a quantum simulator, which by definition should be 
an extremely clean and controllable setup, one can exploit the 
quantum parallelism by an auxiliary degree of freedom [51]. 

In principle, some randomness could be introduced by ran- 
domly varying the lattice constant and the trapping frequen- 
cies within the microtrap array. In the limit of large arrays, 
these would lead to an off-diagonal bond disorder and diag- 
onal site disorder, respectively. In this subsection, we elabo- 
rate on two alternative directions to widen the applicability of 
the phonon-based quantum simulator by introducing disorder 
regardless of the size of the microtrap array. In both cases, 
we shall make use of the electronic energy levels of the ion, 
{| ti)^ I ii)}' to introduce randomness in the phonon Hamilto- 
nian and mimic the effects of disorder. We emphasize that this 
setup shall allow us to control the two usual types of disorder 
independently, namely, diagonal and off-diagonal disorder. 



1. Off-diagonal bond disorder 

We discuss how to induce randomness on the phonon tun- 
neling by a slight modification of Eq. (35), so that the peri- 
odic driving will be responsible for both the disorder and the 
synthetic gauge fields. In the derivation of the driving, we as- 
sumed that a proper choice of the laser intensities, detunings, 
and polarizations would lead us from Eq. (32) to the desired 
expression (33). This condition can be modified so that the 
periodic driving becomes spin dependent. In fact, by setting 

= {^il^y^^^^/2A = -(^lie^)*^lieV2A, one obtains 

Hl = ^^^LliO-fe^^^'^-^i-^^) +H.C., (45) 

where we have introduced cr5 = | | — | |. This ex- 

pression corresponds to a differential Stark shift between the 
two electronic levels. In the regime where cOi^^cOa, this term 
generalizes Eq. (35) to the following spin-dependent driving 
of the trapping frequencies 

ifL - -^lLl^r7^cos(Ak-r?-(OLOcTf/7^^.Z7^., (46) 

where the same constraints (36) over the system parame- 
ters must be fulfilled. Since no other spin operators are in- 
volved in the Hamiltonian, one can treat as c-numbers 
Oi G { — 1, 1}, and carry out the same analysis made in Sec. II. 
In fact, one should simply modify Hamiltonian (39) to ifgff = 

i:|^|ifeff({ty})|{tT})({a}|, where ifeff({tT}) =/reff({CT}) + 
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Veff , such that 
K.n{{o}) = II/^ij(cTi,aj)e'-^'4'^'-.A.^t .^^ 

a i>j 



-H.c, (47) 



and the tunneling amplitudes account for the different spin 
configurations Oj) = j)(T]d(Ti, r]dOj, A^y). In 

Fig. 7, we represent the dressed tunneling amplitudes between 
nearest-neighbors along the direction of the gradient. It can 
be observed that depending on the spin state, one obtains dif- 
ferent values. In particular, for T]^ ~ 1 and ^ 37r/2, the 
tunneling amplitude can attain four possible values 



ii . 



n , 



^0.5. 



(48) 



If we now consider the following initial state po = |^s) (^s | ^ 
p^i^, where |^s) = 11{g}C{g}\{<^}) and p^^ are arbitrary spin 
and phonon states, its time-evolution is 



Pph(0=tr.{p(0}^ 



{<y} 



-///eff({(T})? +///eff({a}> 



(49) 

Note that due to the superposition principle of quantum me- 
chanics, the phonon state explores simultaneously the differ- 
ent tunneling paths i ^ /(OiOj) ^ j with a probability that de- 
pends on the initial spin configuration = p. In fact, 
the measurement of a phonon observable Oph corresponds di- 
rectly to the statistical average over all such tunneling paths 



(Oph(O) = i:/^{CT}tr{Ophe--^eff({CT}>p0^e+^-//e,({CT}).| 



(50) 



which can be understood as the average over all possible real- 
izations of the bond disorder. Therefore, our phonon-based 
simulator can explore the physics of interacting disordered 
bosons in a lattice pierced by an external magnetic field. 



2. Diagonal site disorder 

We now address a scheme to introduce randomness in the 
on-site energies [22]. The origin of such terms is independent 
of the periodic driving, and requires a pair of additional laser 
beams in the Raman configuration (31). In order to distin- 
guish them from the previous laser beams, we will be denote 
them with an overbar. The main difference with respect to the 
previous case is that the beatnote is tuned close to the transi- 
tion between the internal states (b\^ = a)\ — ayi ^ (Oq. In this 
case, the laser-ion Hamiltonian becomes 



(51) 



where we have introduced cr^ = | ti)(ii I- In analogy to the 
derivation of the periodic driving, we express the position in 
terms of the local phonon operators, and expand the Hamilto- 
nian for a small Lamb-Dicke parameter. If the laser beatnote 
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Figure 7: Spin-dependent modulation of the tunneling: Con- 
tour plot of the real and imaginary parts of the modulation amplitude 
^r{ii-ji) the tunneling between nearest-neighbors i ^ j, such 
that ii = 71 + 1, as a function of the driving parameters T7d,A0ij for 
r = \. (a,b) Tunneling in the configuration of anti-parallel spins 
^ii^j ^ {ti^it}^ (c^d) Tunneling in the configuration of parallel 
spins Oi, Oj G {tt?ii}- Also shown in black- white circles the region 
of interest for T]^ ^ 1, A0ij = 37r/2. 



is tuned as follows CO]^ ? 
red- sideband excitation 



COq — COa, one obtains the so-called 



^^^■'io.-^b .e- 
1 a.i 



H.C, 



(52) 



where the bare detuning is = (Oq — (Oa — ^ (Oa^ and 
we assume that <C (Oa in order to neglect the remaining 
terms of the Taylor expansion. In the regime where the laser 
beams are weak enough ^^L^a 5l, it is possible to find the 
following laser-ion Hamiltonian in perturbation theory 



45l ■ 



(53) 



This term arises due to a second-order process where a phonon 
is virtually excited and reabsorbed by a single ion, and leads 
to a differential Stark shift of the atomic levels that depends 
on the number of phonons. When incorporated to the effec- 
tive description of the PAT Hamiltonian (39), it modifies the 
Kinetic energy term. If one is interested in the time span 
t ^ l/^c;ij. th^ dynamics of the spins can be safely ignored, 
and (7^ can be treated once more as c-numbers Oi G { — 1,1}. 
We note that the typical time-scales for the spin flip-flop dy- 



namics would be 



s;ij 



d 0-1 00 Hz, whereas thevibrational 



couplings lie in the /c;ij ^1-10 kHz. Accordingly, the Ki- 
netic energy term of Eq. (39) must be modified to //gff ({^}) = 
^eff({^y})+Veff, where 



a i>j 



a,i a,j 



^ .b .- 



H.c. 



(54) 



In this situation, it is not the tunneling strengths which depend 
upon the spin state (47), but rather the on-site energies. No- 
tice that the strength of these on-site energies must be much 
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smaller that the trapping gradient £a <C AcOa in order not to 
affect the PAT scheme. According to Eq. (50), the system ex- 
plores simultaneously all possible values of the on-site ener- 
gies with probabilities that depend on the spin configurations 
of the initial state. Therefore, the measurement of phonon 
observables yields directly the statistical average over all pos- 
sible realization of £j{{G}) G with a probability 
distribution /?|^| = 

Let us now comment on the extended possibilities of our 
quantum simulator due to the engineered disorder in Eqs. (47) 
and (54). The diagonal site disorder leads to the well-known 
Anderson localization in the non-interacting limit [53]. This 
phenomenon is due to the interference of the different paths 
associated to the scattering of the particles from the random 
on-site fluctuations, and gives rise to exponentially-localized 
wavefunctions and absence of diffusion. The combination 
of Anderson Localization with strong interactions, which is 
also a well-studied problem [42], leads to interesting insulat- 
ing, yet gapless, phases such as the Bose glass. In the case 
of strong bond disorder, a different gapless insulator known 
as the Mott glass arises, which consists of disconnected su- 
perfluid regions of random size [52]. Besides, our quantum 
simulator has the potential of combining both bond and site 
disorder, and tuning them independently, which may pave the 
way towards other exotic insulating phases. In addition to the 
aforementioned Bose and Mott glasses, the simulator can ex- 
plore the random- singlet glass where the bosons form delo- 
calized random pairs [54]. The possibility to explore higher- 
dimensions, longer-range tunnelings, and the effect of syn- 
thetic gauge fields makes our scheme a very versatile tool. 



F. Decorated synthetic gauge flux lattices 

In this section, we describe an additional feature of our 
quantum simulator: the possibility to decorate the lattice with 
any desired pattern of synthetic fluxes. It is thus possible to 
engineer highly inhomogeneous synthetic gauge fields, even 
reaching inhomogeneities at the unit-cell limit. The idea is 
to use the spins to decorate the array with different fluxes by 
exploiting the differential phonon-dependent Stark shift (53). 

We consider a situation where the gradient of the trap- 
ping frequencies vanishes AcOa = 0, so that the regime is dif- 
ferent from that of the site-disorder case considered above 
£a <C AcOa. When the periodic driving frequency (35) and the 
Stark shift (53) fulfill the resonance condition CO]^ = 2eoc/r, 
where r is the integer representing the umber of photons in- 
volved in the PAT, the assisted tunneling will give rise to a 
different phase depending on the spin states of the two neigh- 
boring ions. In order to find the correct expression for the 
tunneling, we readdress the derivation of Sec. II for this par- 
ticular spin-dependent situation. Since we are interested in the 
phonon dynamics, the spins are effectively frozen, and we can 
consider as c-numbers Oi G { — 1,1}. Hence, the dressed 
tunneling (13) must be modified as follows 



where we have introduced /(Oi, Oj) = r(Oi — Oj)/2. By con- 
sidering the possible spin configurations, we find 



(iG\) f ^r(r]d,r]d,A0ij)e-'5(^i+^j), if di = -aj 



^_,(r7d,r7d,A0ij)e+^-5^^i+^j),if(^ 
^o(^d,^d,A0ij), if Oj 



: 1 

-CTi = 1 

ai = ±l. 
(56) 

Interestingly, the phase of the tunneling between sites j i 
depends on their internal spin state Oj ^ (Ti. If their spins are 
parallel Oj = Oj = ±1, the phase vanishes and thus the tun- 
neling does not contribute to the synthetic gauge field. Con- 
versely, when the spins are anti-parallel, the phases contribute 
to the gauge fluxes with a sign that depends on the particu- 
lar spin ordering. We have calculated the consecutive phonon 
tunneling around a square plaquette r? ^ r? ^ ^ rj^ ^ r? 
for all the possible spin configurations 



.(1) 



^^^)\eM{<y}) for 



^d;ij 



(^i+^j), /^ij=/^j^/(.,.,)(T]d,r]d,A(/)ij) 



(55) 



All the possible encircled fluxes 
r = 1 have been represented in Fig. 8(a), where we observe 
that there are 9 possible fluxes out of the 2^ = 16 possi- 
ble spin configurations. These fluxes correspond to G 
{O,±0i,±02,±0+,±0-}, where 0,- = Akidi, and 0± = (0i ± 
02) /2. Accordingly, we have 9 different tiles that can be 
used to decorate the underlying microtrap array with a par- 
ticular distribution of fluxes. Let us emphasize that the partic- 
ular distribution of tiles is completely determined by the spin 
state l^s) = , CTiv), which can be initialized at will in 

trapped-ion experiments. Let us also note that we could also 
exploit this result to introduce randomness in the gauge fields 
by considering a linear superposition of the spin states. 

In Fig. 8(b-c), we represent two-color flux lattices that cor- 
respond to a staggered magnetic field along both principal 
axes. In Fig. 8(d-e), we represent three-color flux lattices 
where a staggered flux alternates with a vanishing flux. In 
Fig. 8(f), we represent a checkerboard flux lattice, and finally 
in Fig. 8(g), we represent a limiting case of a six-color flux lat- 
tice, which is very interesting from a physical point of view. 
By setting 0i = tt, and 02 = ^, the two fluxes ±0i = ±7r = 
7rmod27r are equivalent, and lead to a homogeneous TT-flux 
model in a square lattice. Additionally, the four remaining 
fluxes vanish for this choice ±0+ = ±0_ = 0mod2;r. There- 
fore they contribute with a local defect over the TT-flux lattice, 
which could bind excitations with anyonic statistics when the 
longer range tunnelings are taken into account [55]. 

For each of the decorated flux lattices, it is possible to find a 
particular inhomogeneous synthetic gauge field As(Oi,Oj) so 
that the effective phonon Hamiltonian is rewritten in a stan- 
dard Peierls form 

a i>j 

(58) 

The particular pattern of the spins will determine the inhomo- 
geneous gauge field, and the way the lattice is decorated with 
fluxes. Let us also note that the complexity of Fig. (8) will in- 
crease when the larger plaquettes due to long-range tunnelings 
are also taken into account. 
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Figure 8: Decorated synthetic flux lattices: (a) Unit cell for a square array of microtraps with all the possible configurations of spin states 
G[ G {1,-1}, which are represented by circles at the vertices of the cell. According to the spin-dependent enclosed flux, we have a total of 9 
different plaquettes that can act as tiles to construct the decorated flux lattices, (b-c) Staggered flux lattices, (d-e) Staggered flux lattices with 
alternating strings of zero-flux plaquettes. (f) Tetra-flux checkerboard lattice, (g) TT-flux lattice with a localized defect consisting of a zero-flux 
plaquette. 



IV. SUMMARY AND OUTLOOK 

We have introduced the two key ingredients to realize PAT 
experiments in micro-fabricated ion traps. The first ingre- 
dient is a gradient of the trapping frequencies achieved by 
the local control of the trap electrodes. The second corre- 
sponds to a Raman-beam configuration, which presents dif- 
ferent regimes that provide (i) the periodic driving of the 
trapping frequencies, (ii) the on-site phonon-phonon interac- 
tions, (iii) the bond/site disorder, and also (iv) an exotic flux 
decoration of the microtrap array. We believe that such in- 
gredients are at reach of current microtrap technology, and 
their correct combination will give raise to a very versatile 
quantum simulator for many-body bosonic models. Such 
photon-assisted-tunneling toolbox for quantum simulations 
can be summarized in the following general Hamiltonian 
= Lie} M{a}) + Veff) \{G}){{a} I , where 

a i>j 

+ + (59) 

a i 

i,a,7 

where the particular expression for the dressed tunneling and 
synthetic gauge field will depend on the configuration of the 
frequency gradient and the periodic driving. Let us list the 
possibilities that have been explored in this work: 

i) Dynamic localization: The tunneling amplitude /^y = 
/^^j^o(^d7 ^d, does not depend on the spin state, and the 
synthetic gauge field vanishes Af = 0. This is achieved in the 
regime of vanishing gradient AcOf^ = 0, and setting the beat- 
note of the Raman beams (Ol <C 0)^. By tuning r]d, one can 



find a value where the tunneling strength vanishes and thus 
the particles are dynamically localized, a phenomenon also 
known as coherent destruction of tunneling [25]. 

ii) Synthetic Abelian gauge fields: The spin-independent 
tunneling amplitude is = j)(r7d, T]^, A^jj), where 
/(i,j) = r(/i —71) depends on an integer r. The synthetic 
gauge potential is As(r) = — ^oj^i, and follows from the 
regime with a finite gradient Aco^ ^ /^^j, such that the beat- 
note of the Raman laser beams fulfills the resonance condition 
(Ol = AcOfx / r. In this case, phonons behave as charged particles 
that move in a two-dimensional plane pierced by an orthogo- 
nal magnetic field whose flux 02 = A^2<^2 can be modified by 
varying the Raman wavevector. Phenomena typical of integer 
quantum Hall samples [39, 40], or bosonic flux ladders [41], 
can also be observed in this platform. Besides, in combination 
with strong phonon-phonon interactions, one can find bosonic 
versions of the fractional quantum Hall states [43]. 

iii) Synthetic non-Abelian gauge fields: The above scheme 
can be generalized to the non-Abelian gauge group SU(2), 
such that both in-plane vibrational modes play the role of a 
flavor component. We have described in detail a particular 
SU(2) gauge field, which requires the vibrations along each 
direction to be subjected to an opposite frequency gradient. 
Using the same assumptions as in the Abelian case, the tun- 
neling amplitude becomes /^y = ^tUj^/a(iJ)(^da, ^da, A^y), 
with j) = ra{i\ — j\) such that = —Vy = r. Hence, the 
synthetic gauge field A^^ = —B^yTz^i becomes a SU(2) op- 
erator acting in the flavor space. This scheme opens a route 
towards a bosonic counterpart of the quantum spin Hall ef- 
fect [49]. 

iv) Bond and site disorder: The dressed-tunneling ampli- 
tude /^jj ( (Ti , Oj ) = /^^j ^^(i j) ( r]d cJi , r]d Oj , A0ij ) , and the on-site 
energies ea<^\^ take on different values depending on the spin 
configuration. If the initial state is a linear superposition of 
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different spin configurations, the phonon dynamics is deter- 
mined by a random Hamiltonian with bond and site disorder. 
This regime is achieved for a gradient AcOa ^ J^^ , and laser 
beatnote (Ol = AcOa/r giving rise to a state-dependent periodic 
driving. Besides, an additional Raman beatnote tuned close to 
the atomic transition G)[^^ (Oq — COa gives raise to a phonon- 
dependent Stark shift in the limit of large detuning. In the non- 
interacting regime, this tool allows us to explore the physics 
of Anderson localization [53]. By adding strong interactions, 
it yields gapless insulating phases such as the Bose glass [42], 
the Mott glass [52], and the random- singlet glass [54]. 

vi) Decorated flux lattices: In the absence of the fre- 
quency gradient, one can tune the Raman lasers beatnote in 
resonance to the above phonon-dependent Stark shift (Ol = 
lEa I Once again, the dressed- tunneling amplitude becomes 



spin-dependent J^.. (di, Oj) = J^i^^f{Gi,Gj){Vd, ^d, A(/)ij) where 
/(Oi, (Jj) = r((7i — Oj)/2. Moreover, the synthetic gauge field 
also depends on the spin configuration and we can make dec- 
orated flux lattices as those shown in Fig. 8 by selecting a 
particular spin state. Phenomena related to charged particles 
under inhomogeneous magnetic fields can be explored, such 
as staggered fields or ;r-flux lattices with defects. 
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